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THERMAL DYNAMIC MODELING STUDY 
by I. U. Ojalvo 

Grumman Aerospace Corporation 
SECTION 1 INTRODUCTION AND SUMMARY 

The space shuttle thermal design must accommodate the environments associated 
with prelaunch cooling, boosted flight, earth orbit operation, reentry, and postlanding 
soak-back. The design conditions which are believed to be of importance to the 
shuttle's thermal-structural response are: 

• Reentry - which will establish the most severe thermal environment, thus 
designing most of the thermal protection system (TPS) and causing the 
greatest primary -structure thermal stresses 

• Abort maneuvers - during flight these are of less importance because of the 
present decision to define an "admissible" abort as one which a reentry 
based design could survive. Pad aborts, however, in which the orbiter 
accelerates away from a booster which is about to explode, are concerned 
with tradeoffs between high separation acceleration loads, and thermal- 
dynamic loads resulting from such an explosion 

• Stage separation - engine firing of one vehicle in close proximity to portions 
of another vehicle's structure, such as a fin, could cause high local thermal 
stresses. This, of course, is highly dependent upon the structural con- 
figuration of the mated vehicles and fuel tanks 

• Orbital environment - the vehicle temperature prior to reentry strongly 
governs its internal temperature during reentry. It is possible that different 
relative solar positions, prior to reentry, could account for significantly 
different thermal states 

Requirements associated with the above phases present a wide range of formidable 
thermal-dynamic structural problems. Although related analysis procedures have improved 
significantly in the past decade, analytical difficulties still exist due to geometric 
nonlinearities, inelastic and nonlinear material behavior, complex liquid-solid 
interactions, and the presence of poorly defined damping mechanisms. Thus, 
analytical solutions must be supported by carefully controlled laboratory tests. 
Furthermore, the large size of the vehicles involved requires that much of the 


experimental work be performed on reduced-scale models. However, before mean- 
ingful tests can be conducted, it is necessary to understand the proper thermal dynamic 
similarity laws and to evaluate the feasibility of their implementation in a practical 
modeling program. The present work is directed towards providing basic information 
for designing the models and conducting the required tests. Items considered are the 
development and interpretation of thermal dynamic structural scaling laws, the iden- 
tification of major related problem areas, and the presentation of viable model fabrica- 
tion, instrumentation, and test procedures. 

A device employed to aid this general stucfy was the selection of a particular 
component — an orbiter fin — to serve as a focal point for portions of this effort. 
Thermal structural and stiffness analyses were performed upon the fin for a typical 
reentry maneuver and for a l/10th reduced scale replica model of it. In addition, 
manufacture, test, and cost analyses were also performed upon the fin model. 

Thermal and structural dynamic similarity laws are derived and interpreted 
in Section 2 with regard to general shuttle design concepts. It is shown that the 

l 

major difficulty associated with application of the scaling laws is that the internal 
radiation, convection and conduction modes of heat transfer do not scale uniformly 
on a reduced scale replica model. For example, the ratio of model to prototype 
time-scale is proportional to length-ratio squared for conduction, but length-ratio to 
the first power for radiation and convection. A second difficulty is that, compared to 
the prototype, the model time-scale may shrink so substantially as to preclude the 
use of conventional steady-state vibration testing. 

Model design of the orbiter fin selected for detail consideration is reported in 
Section 3. Geometric distortions, consistent with practical manufacturing require- 
ments, are discussed. Results of the study indicate the feasibility of replica model 
manufacture but underscore its extremely high cost. 

Temperature, stress and stiffness analyses for the fin and its model undergoing 
a reentry maneuver are presented in Section 4. Geometric distortion and replica size- 
scaling errors in the model are quantitatively evaluated and found to be small. 
Furthermore, it is demonstrated that thermal stresses may be greatly relieved 
through implementation of intelligent design concepts. In particular, for the case 
studied it is shown that the use of corrugated (as opposed to noncorrugated) spar 
webs reduces primary thermal stresses dramatically. Thus, reductions in 
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cross-sectional twisting stiffnesses (up to 8%) were not influenced by destabilizing 
thermal stresses, but were due almost entirely to reductions in material stiffness 
properties associated with the elevated temperatures. 

Dynamic testing at elevated transient temperatures is discussed in Section 5. 
Such items as experimental methods, instrumentation, excitation methods available, 
and necessary data processing techniques are identified. 

Concluding remarks are presented in Section 6. Recommendations are offered 
regarding the need for conducting reduced-scale dynamic tests of major structural 
components at elevated temperatures, including practical suggestions for implementin 
such tests. 

It is hoped that the present "paper-study" will alert the aerospace community to 
some important aspects associated, with not only shuttle -related models, but with 
heated structural models in general. Another goal is that this effort provide useful 
information towards planning meaningful thermal dynamic experiments. It is strongly 
felt, however, that only through careful experiments, with actual full-scale and 
model hardware, will a true appreciation of the need for such tests become evident. 
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SECTION 2 SCALING LAWS 


A number of works on general modeling laws, structural dynamic similitude, 
or thermal scaling alone have been presented, several of which are included in 
References 1 through 3. The present work differs from these, however, in that an 
attempt is made to identify and interpret the relevant thermal dynamic scaling laws as 
they apply to space shuttle primary structure. Also included is an investigation into 
the importance and feasibility of their application in a practical modeling program. 

The basic dimensionless groups for thermal dynamic model studies of aircraft 
and spacecraft structures are derived from equations governing internal heat transfer 
and structural dynamic response. Although other derivation procedures are available, 
this method is preferred since it provides more insight than the mere knowledge of 
the dimensionless variables entering the problem. 

After deriving rather general scaling laws between model and prototype, they 
are specialized for application to space shuttle primary load-carrying structure. 
Parameters associated with the possible use of different temperatures, materials, 
and method of construction between the prototype and the model are also discussed 
in light of an approximately 1/1 Oth reduced-scale model. 

Although the aerodynamic environment of the shuttle prototype is the primary 
source of heating, the possibility of testing a 1/lOth scale model in a high stagnation 
temperature wind tunnel is not likely because of the great expense required to develop 
such a facility. In connection with this, it is noted that the similarity laws for com- 
plete aerothermoelastic testing were investigated by Dugundji and Colligeros (Re- 
ference 4), who found that similitude in the general case was difficult to achieve for 
a scale ratio other than unity. Their conclusions suggest the use of "incomplete 
aerothermoelastic testing" and "restricted purpose" models. For practical reasons 
it is anticipated that model heating would be best accomplished by the use of quartz 
lamps programmed to generate predetermined temperature-time profiles at selected 
points instrumented with thermocouples. Also, since the thermal protection system 
(TPS) will not be part of the load carrying structure, only TPS backface temperatures 
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and TPS mass need be simulated, thus simplifying model fabrication and lowering test 
temperature levels. 

2. 1 THERMAL SCALING LAWS 

2.1.1 Derivation 

The mechanisms of internal heat transfer which may be present simultaneously 
are: conduction, convection, and radiation. Combining these into a single term (q), 
defined as the rate of heat absorbed by a body, and applying the first law of thermo- 
dynamics (conservation of energy), yields: 

'C V = q (1) 

where P, c, V and T, respectively, are the density, specific heat, volume, and tem- 
perature of the body, and t is time. 


Applying Eq. 1 to a given element i of a lumped thermal node idealization, and 
expanding the right-hand side explicitly into its various components , results in the 
internal heat balance equation: 


dT. ^ A. r-' 

P. c.V. 1 = > k. -— iHL(T - T. ) + / h. A. (T - T.) 

i i i im LT ' m r L in in ' n l 

° m im n 


(conduction) 


(convection) 


+ / cr A. F. 
£-> lr lr 
r 


— 4 




+ ^ q . A. 

s is is 


(radiation) 


(heat flux source) 


(2) 


The ! specific form used in Eq. 2 is not unique, but represents the formulation 
typically programmed and employed for complex structural temperature analyses. 
The variables involved are defined as follows: 

k. , L. , A. are the conductivity coefficient, length and area, 

im im im ’ 

respectively, of the heat conduction path joining 
elements i and m 

h. , A. are the convective film coefficient and the surface area, 

in in ’ 

respectively, between elements i and n 
<r is the Stephan-Boltzmann constant 
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A. is the surface area of element i that "sees" element r 

lr 

F. is a radiation heat transfer coefficient which depends on the 

geometry and surface emissivities of the structure 

q. is an additional heat flux rate associated with boundary- 

condition nodes 


A. is the nodal surface area through which additional heat 

is 

flux passes to element i 


T., T , T , T are absolute temperatures 

1 m* o r 

The temperature and thermal time scales, 0 and T T » 

by the ratios 

T' t * 

6 =~ and T = 7 - 
X It 


respectively, are defined 

(3) 


where the primed quantities refer to model quantities*, and the length, area and volume 
scales are defined as follows: 



Subscripts are used on each A to accommodate the possibility of employing distorted, 
as opposed to geometrically replica, models.- For the case of replica model scaling, it 
is obvious that the L, A, and V subscripted length-ratios are all equal. 


Applying Eq. 2 to both prototype and model structure, employing Eq. 3 and 4, and 
comparing the results yields the following scaling laws to maintain the same scaled 
temperature-time profile in both: 2 

k k’ T 


(conduction) 


Ac 


pc p'c' 


X L X V 


(5) 


(convection) 


_h_ 
p c 


h’ 


p'c' 


T T X As 


(Cont. on following page) 


*A distinction is made between thermal time scale ( Tp) and structural dynamic time scale 
.(to be introduced in section 2 . 2 as rp) since thermal events occur at a much slower 
rate than vibratory mechanical events for the problems under consideration. Thus, for 
engineering analysis purposes, thermal states may be assumed quasi-static with regard 
to dynamic response, which for modeling purposes implies that T-p may be selected in- 
dependent Of Tq. 
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(radiation) 


(heat flux) 


F F - r T^ X As 

pc P'c' X 3 y 

2 

J*_ = T T X As 
P C P ' c* x 3 v 0 


(5) 

(Cont. from previous page) 


where primes over quantities refer to model parameters, and indices im, ir and is have 

2 

been eliminated for simplicity. However, subscripts c and s have been added to the \ A 
to denote scaled area ratios associated with heat transfer through either a cross-sectional 
or surface area, respectively. Combining Eq 5, and using Eq. 3 and 4 yields the dimension- 
less groups 


kt 


ht 


crFtT 3 


and 


qt 


(6) 


PcLr P cL ’ pcL P cLT 

where the Stefan-Boltzman constant, <x, and the absolute temperature, T, have been in- 
troduced to make the last two groups dimensionless, and L is a characteristic length (not 
necessarily aphysical distance between two points). The first of these dimensionless groups 
is the Fourier Number, while the second and third are products of the Fourier Number with 
the Biot and Stefan numbers, respectively (Reference 5). 

2.1.2 Interpretation 


It can easily be seen that if the model is a geometrical replica, is made of the same 
material, has the same surface finish and heat flux rate as the prototype, and internal 
convection coefficients are identical, it will not be possible to satisfy all the require- 
ments for any length-scale other than unity since Eq.5 becomes: 

2 

(conduction) r^, = A 


(convection) r^, = X 

(radiation) = A/0 ^ 


(heat flux) = X 0 


whose only nontrivial solution* is 


(7) 


X = 


= e = i 


(8) 


*The trivial solution referred to is X = r T =0, which is not physically useful. 
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If conductive heat transfer does not enter into the problem (or can be disregarded 
because it is small) the first of Eq. 7 need not be satisfied and the remaining two admit 
the solution 

X = r T and 0 = 1 (9) 

However, if conduction is the only mode of heat transfer to be considered, then only the 
first of Eq.7 need be satisfied. 

It should be noted that, the specific heat, c, and density, p % are material pro- 
perties while the convective film coefficient, h, and the radiation coefficient, F, depend 
upon other characteristics such as surface roughness, and surface coating. Although the 
conductivity is also a material property, the conductivity coefficient, k, may refer to 
an overall conductive path and include such considerations as joint conductivity, which 
is obviously a function of the type of joint. 

It may then be possible, even when using geometric scaling and the same material 
and temperatures for model and prototype, to satisfy all the requirements of Eq.5 (if 
not exactly, then with sufficient approximation) by making the model in such a way that: 



Jv_ _ q' X (10) 

h F q r T 

It should be noted that in all probability there will be difficulty in geometrically 
scaling all member dimensions accurately (such as the cross-section of a thin post). In 
such cases, reference must be made to Eq. 4 and 5 to identify the pertinent physical dimen- 
sionless ratios so as to evaluate any distortions introduced by the model. 

The basic conflict in accurately modeling all internal modes of heat transfer, as 
applied to shuttle primary structure, may be investigated in a preliminary way as 
follows. For discussion purposes only, consider the problem of modeling a wing with 
approximately 0.9 m (3 ft) between upper and lower titanium covers, 315°C (600°F) 
bottom-side and 204°C (400° F) top-side temperatures and negligible chord and spanwise 
thermal gradients. 
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The black-body radiation flux corresponding to 315 C and 204 C (600°F and 400°F) 

is about 10 and 4 — ~ ca * (2200 and 940 Btu/hr-ft^), respectively. If hot and cold 

cm 2 -min g-cal 

regions can "view” each other, net radiant flux levels will be of the order of 6 — % . 

cm‘- min. 

(10^Btu/hr-ft2). By comparison, the conduction heat transfer for a gradient of lllc° 

over a 0. 9m (200F°over 3 ft) distance for titaniumat a wing area density of 1 gm/cm (21b/ 

g-cal „2. 


ft^) would only be on the order of 0. 006 - 


(1 Btu/hr-ft ) (it can also be shown 


cm^-min 

that internal natural convection effects for this case are also small). Hence, the 
principal influence on the structural temperature distribution would be the net radiant 
flux. Moreover, since the structural elements, exclusive of any honeycomb construc- 
tion, have relatively thin gages (say 0.6 cm (1/4 in. ) maximum for the prototype), it 
can be shown that the temperature variation within an element will be small across its 
thickness. Therefore, it appears reasonable to neglect transient conduction considera- 
tions in a first approximation to thermal modeling and select r = X as the appropriate 
time scale. 


A similar result was obtained for an ablator or reusable surface insulation (RSI) 
protected aluminum wing using a 48 node thermal model (see Fig. 1) with ablator 
backfacc heating rates shown in Fig. 2. Runs were made for the prototype and a 
l/l0th scale replica model. Results for nodes 1 and 11 (Fig. 1) are plotted in Fig. 3. 
They reveal peak temperature differences of 85C° (150F°) and maximum thermal 
scaling errors of only 8.5C° (15 F°). Greater scaling errors occur, however, for 
nodes with larger conduction couplings. The largest discrepancy encountered with the 
present idealization was 28C° (50F°), which occurred at node 35 (Fig. 1 and 4). It is 
felt that the errors indicated for this node are overly conservative (i. e. , they appear 
higher than they really would be in an actual test) as the mass concentration used was 
high and all radiation couplings to this node were neglected. Nevertheless, the over- 
all results are consistent with what might be expected when dealing with lower struc- 
tural temperatures and higher material conductivities. Both of these factors tend to 
increase the ratio of conduction to radiation heat transfer, thus causing greater, but 
acceptable, modeling distortions. 

A more detailed thermal analysis upon a high temperature alloy structure with- 
out a TPS is presented in Section 4.1. As might be anticipated, the percentage errors 
associated with it are lower because the ratio of radiation to conduction is even higher 
than the two cases discussed previously. 
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Fig. 3 Comparison of Shuttle Wing-Segment Skin Temperatures for Aluminum Prototype 
and 1/10th Scale Aluminum Model 
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Fig. 4 Comparison of Shuttle Wing Hardware-Attachment-Point Temperatures 
for Prototype and 1/10th Scale Model 
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2.2 STRUCTURAL DYNAMIC SCALING LAWS 
2.2.1 Derivation 

Elasticity equations governing the stresses and large deflections in a three- 
dimensional, isotropic, heated body which obeys Hook’s law can be written in indicial 
tensor form as 

<>\ 

■* ' '*> s- 1 '' 1 - 2 ’ 3 (ii) 

^-5., cr +a« T (12) 
E jk rr ]k v ' 

where 


and 


/du au k 

e = 1/2 + — 


du du 


jk 


v S x 


k 


ax. 


dx. 


d \ 


r) - 


1 +v 


jk 


cr , - cr, 
Jk kj 


€ • 


jk " *kj 




u. 

J 


t 

g j 

V 

E 


*jk 


are the stresses 

are the strains 

are the cartesian coordinates 

are the displacements 

is the mass density 

is time 

are components of the gravity or body-force field 

is Poisson's ratio 

is the modulus of elasticity 

is the coefficient of thermal expansion 

is the temperature above the stress-free state 

is the kroneker delta, which is unity when the subscripts 
are equal and zero otherwise 


In addition to the three dynamic equilibrium and six symmetric stress deflection 
relationships (Eq. 11 and 12, respectively) in nine unknowns (three deflections and six 
stresses), it is necessary that stresses and/or deflections be specified on the bounda- 
ries in an appropriate manner. 
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Length, time*, and temperature ratios are defined as before and a deflection 
scale, A u , is added to this group to give 


X 

x 


x' 

X 



- T' 

8 = -=r and X 
T u 


u / 

u 


(13) 


where primes again denote model parameters. 

Applying Eq. 11 and 12 to the model and substituting Eq. 13, lead to the follow- 
ing requirements for similarity: 


(dynamic equilibrium) — = ~ — 9 — 

P 


and 


(stress-deflection-temperature) 


•= 1 , 


v' 


X 


' E »v 


, o’ cr' . 

= 1 , = , « = «0 

E E' 


(14) 


Combining Eq. 13 and 14 yields the six dimensionless groups 


pLu pgL 

9 * 

at & 


U y <T 

IT* » T* and aT 


(15) 


It should be noted that the material properties, E, a , p and v may be functions 
of temperature, T, which may itself vary with time. This thermal variation is generally 
much slower than the variation of deflection with time, thus justifying the distinction used 
between "thermal" and "dynamic" time-scales, and r^, respectively. 


Because the analysis of aerospace structures may be simplified by considering 
them to be composed of bars, beams, membranes and plates, it is useful to consider 
the thermo-structural-dynamic equations of these elements as well. As an illustra- 
tion of the derivation techniques involved, the beam element is considered in detail 
in the Appendix. However, the appropriate scaling parameters for several other 
elements are only summarized in Table 1. 

2.2.2 Interpretation 

The first of the preceding dimensionless groups relates to dynamic loads, the 
second to stresses produced by gravity effects, and the third arises in large deflec- 
tion problems. The interaction among these three effects occurs when a preload 

*In this case the time referred to is "dynamic", as opposed to the "thermal" time of 
Section 2.1.1. 
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TABLE 1 SCALING LAWS FOR ELASTIC STRUCTURAL DYNAMIC MODELS 


Element 

Scaled Quantity 

Scaling Laws 


Dimensions and Coordinates 

L' x ' 

L " x ~ A 


Deflections 

— - XL - w ' - \ 
u - v w 


Strain (linear and nonlinear) 


General 

Time 

V _ Ip 7e’ 

t T D V p / E 


Thermal Strain 

g’T 1 = 
a x ~ 


Frequency 

V ’ 1 /E’/p' 

f \ y e/p 

Bar & Beam* 

Stretching Stiffness 
Mass Per Unit of Length 

(EA) / _ _E^ x 2 
(EA) E 

(PA)' _ P' x 2 

( *>A) p 


Bending Stiffness 

(El)' E' 
(El) ~ E 


Twisting Stiffness 

(GJ)' _ G'^4 
(GJ) G 


Stretching Stiffnesses 

(Eh)' _ E' 
(Eh) ~E 

Membrane & 
Plate* 

Shearing Stiffness 

(Gh)' G' 
(Gh) ~ G 


Bending and Twisting Stiffnesses 

~d! d' e' 

1 . xy . x 3 


(u' = v) i=x,y 

D. D E 

i xy 


Mass Per Unit of Area 

(Ph)'_ P\ 
(Ph) P 

♦Neglecting shearing deformation 
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stress state, caused by gravity, influences the dynamic characteristics through a non- 
linear stiffness interaction. One way of seeing this directly is to consider a replica 
model and prototype made of the same materials. Under these assumptions Eq. 14 
permits the solution: 



x = g/g' 


(16) 


The significance of Eq. 16 is that the model and prototype strains and stresses 
are equal if the length scale equals the time scale as well as the ratio of gravitational 
fields between prototype and model. This latter condition on gravity field ratio is 
highly restrictive in selecting a model scale, but its neglect can cause significant 
modeling errors. Just such an effect was noted by Grimes, et al (Reference 6), 
who reported that joint flexibilities of a l/lOth scale model of the Saturn V were 
seriously affected because gravitational preload effects were not modeled properly. 

Another important type of preload can occur with thermally induced stresses 
which are associated with the dimensionless groups <*T and a /E. As may be seen 
by Eq 16 (ignoring the last one on g/g' ), such effects are accurately simulated in a 
replica model if the temperatures between model and prototype are identical. 


2.3 MODELING DISTORTIONS 
2.3.1 Material Considerations 

Parameters associated with the use of different temperatures and materials have 
been investigated. Observation of the dimensionless groups associated with internal 
heat transmission, Eq. 6, and thermo-structural-dynamics, Eq. 15, reveals all the 
pertinent material properties involved in a related modeling study. 

For reference purposes, Table 2 reviews the three popular orbiter design con- 
cepts, the Aluminum-Ablator or RSI/heat-sink, Titanium/TPS-protected, and Nickel- 
base superalloy /hot-structure, their temperature design limits, and dominant modes 
of internal heat transfer. The aluminum concept presents the biggest thermal modeling 
problems because multiple significant modes of heat transfer with conflicting scaling 
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laws are present. However, heat conduction in primary structure is still considered to be 
secondary to radiation energy transfer even for a heat-sink design. Furthermore, the 
al uminu m concept is associated with the lowest thermal gradients and so involves the small- 
est thermal-dynamic interaction. For these reasons, and also because radiation effects can 
be modified somewhat through application of surface coatings to change emissivity, it ap- 
pears that density and specific heat are the major thermal parameters to consider in select- 
ing a model material. 


TABLE 2 SHUTTLE ORBITER DESIGN CONCEPTS 


Orbiter Design 

Primary Structure 

Thermal 

Dominant Modes of 

Concept 


Design Temperature 

Protection 

Heat Transmission 

Heat Sink 

Aluminum 

175°C (350°F) 

Ablator 
or RSI 

Radiation & Con- 
duction 

TPS Panels 

Titanium 

315°C (600°F) 

Metallic & 
Nonmetalic 

Radiation 

Hot Structure 

Inconel & 
Rene 

540 (1000) & 
650°C (1200°F) 

None 

Radiation 


The key variable to consider in shuttle structural dynamic simulation is the modulus 
of elasticity, even though the coefficient of thermal expansion affects stiffness (through 
prestress level) and density influences a system's dynamic response. The reason for this 
is that the basic stiffness parameter, E, is temperature dependent (see Fig. 5) and temp- 
erature varies with time throughout the structure. This transient effect is accentuated in 
a 1/lOth scale model since the time scale also shrinks to 1/lOth if conduction effects are 
neglected. Therefore, in a transiently heated structural dynamic model it appears desir- 
able to employ the same material and temperature levels for the model as for the prototype, 
particularly since the primary structure temperature levels are achievable in the laboratory 
(see Section 5.2); and dynamic instrumentation, which function acceptably at these thermal 
levels, is available (see section 5.2). 

2.3.2 Geometric Considerations 

The degree of geometric similarity which is maintained between any model and its 
prototype will influence the applicability, quantity, and type of information which may be 
obtained from test. Roughly speaking, the need for geometric similarity in a model is 
inversly proportional to the understanding of the physical phenomena involved. Thus only 
when the effect under investigation is fairly well understood can distorted models be useful. 
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E q = ELASTIC MODULUS AT ROOM TEMPERATURE 
MATERIAL E n 




Fig. 5 Elastic Modulus Ratio, As a Function of Temperature, for Several Space 
Shuttle Candidate Materials 
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The major difficulties with replica models are associated with scaling com- 
plex joints, attachment fittings, thin skins and mechanical fasteners without ser- 
iously magnifying tolerance-error percentages. Furthermore, the costs associated 
with replica model manufacture, as will be shown in Section 3. 3, can become quite 
significant. 

Therefore, the use of "adequate" models, where possible, is suggested. Illus- 
trations of adequate modeling procedures include: replacing a corrugated spar web 
(acting purely as a shear panel) with a panel which matches its average shear stiffness, 
and substituting a stiffened cover skin with one having a different shape and number of 
stiffeners from the prototype, but possessing properly scaled stretching stiffnesses. 
However, caution must be taken to see that phenomena involving parameters which 
have not been scaled do not occur in either model or prototype; e. g. , improper model- 
ing in the above examples would exist if bending stiffnesses are not scaled and 
buckling occured in either model or prototype. 
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SECTION 3 MODEL DESIGN AND FABRICATION 


3. 1 COMPONENT DEFINITION 

To better evaluate the feasibility of designing, manufacturing and testing a 
thermal dynamic model, as well as obtaining some idea of the cost of performing re- 
lated tests, the vertical fin of an early Grumman-conceived orbiter (designated as the 
H3T) was selected as a focal point for detailed study. The effort described herein is 
a "paper study, " in that an actual physical model was not built and tested. However, 
sufficiently realistic design and analyses have been performed and the necessary 
manufacturing and scaling approximations were evaluated. Thermal dynamic test 
and data-reduction cost estimates to conduct such a program are reported in Section 
5. 2. 5. 

Originally, the component investigated was conceptually designed as both a 
TPS-protected/titanium and an unprotected "hot"/inconel-rene structure. Although 
indications were that the fin with a TPS would be lighter, the unprotected fin was 
selected for detailed analysis since it posed more severe heated primary-structure 
design and testing problems. 

To lend some perspective to the example selected. Fig. 6 presents the overall 
H3T orbiter. A more detailed design of the fin itself is shown in Fig. 7. Basically, 
the primary structure of the hot fin is comprised of a two-cell swept box beam with 
single- sheet/corrugated skin and truss-work ribs which attach to three fuselage ring 
frames at the spar roots. Figure 8 shows the detailed finite element idealization 
used in optimizing the fin for minimum weight and analyzing for thermal stresses. 

6 

The net airload used to design the prototype was a 9. 7 x 10 newton (218, 000 lb) 
transverse resultant, limit load. The assumed pressure distributions, which vary linerary 
in all directions with planform length, are shown in Fig. 9a for the tip and root section. 

The location of the four corners of the resulting center-of-pressure envelope are shown in 
Fig. 9b and, as the pressures are fully reversible, they provide eight different load conditions. 

The projected surface area surrounding each note was computed and a load equal to 
this projected area multiplied by the average pressure at that node was calculated and placed 
at that node. Aerodynamic drag-loads were neglected for this phase of the analysis. As a 
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Fig. 6 Preliminary Design of Grumman H3T Orbiter 
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Fig. 7 Preliminary Design of H3T Orbiter Fin 
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LEFT COVER 


RIGHT COVER 


Fig. 8 Finite-Element Idealization of H3T Orbiter Fin Structure 
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precaution, the rudder was designed fail-safe, i.e. , it consists of two segments, each sup- 
ported by two hinges and an actuator. 

The Grumman Automated Structural Optimization Computer Program (ASOP) (Ref- 
erence 7) was used to automatically size the skins, webs and rib-stud members so as to 
minimize structural weight while not overstressing any member under the above specified 
airloads. Program input to ASOP involved structural idealization, material and buckling 
section allowables, loading conditions and practical minimum gage sizes. 

3.2 MODEL DESIGN 

The philosophy adopted in designing the H3T hot-fin model, similar to that used 
in the l/lOth scale model of the Saturn V (Reference 8), was to scale all dimensions 
by l/lOth and to use the same materials as the prototype where feasibility and reason 
permitted. The resulting design is shown in Fig. 10. Because of certain practical 
limitations associated with the miniaturization process, special manufacturing con- 
siderations and model distortions were necessary. The major ones are reviewed 
here. 

• Rivets - Spot or scam welds in the model replace all rivets used in joining 
the corrugations to the skin face- sheet, rib, and spar-caps 

• Corrugations - The shape of the skin corrugations had to be altered to pro- 
vide adequate size flats to permit welding to the skin face-sheet. This 
caused a significant geometric distortion in the corrugated portions of the 
model. To validate the use of the proposed distortions, the stiffness 

scaling laws required that the depth, pitch, and cross-sectional area-per- 

2 

pitch be closely maintained in the ratios of X , X , and X > respectively, 

(where X = 1/10 is the ratio of model length to corresponding prototype 
length). Figure 11 compares the prototype corrugation's cross-sectional 
shape with that of the model. A separate investigation was performed to 
estimate the skin thermal error introduced by this geometric distortion, the 
results of which are presented in Section 4.1.2 

• Truss Members - The rib struts of the prototype are cylindrical tubes, the 
smallest of which have a 2. 54 cm (1-in.) OD and 0. 04 cm (0. 016-in. ) wall 
thickness. A l/lOth geometric scaling in such cases would be extremely 
costly to manufacture and a revised cross section (shown in Fig. 12) is proposed 
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PROTOTYPE - FULL SCALE MODEL - 10 x FULL SCALE 


Fig. 11 Skin Design Details for H3T Orbiter Hot Fin — Prototype and Model 




FULL-SCALE 

PROTOTYPE 


10 x FULL-SCALE 
MODEL 


Fig. 12 Cross-Section for Rib Truss Braces of H3T Orbiter Hot Fin — Prototype and Model 


29 


• Skin Gages and Tolerances - In attempting to scale the variable exterior 
skin thicknesses, there was no point in attempting to grade them in incre- 
ments of less than 0. 00254 cm (0. 001-in. ). The reason for this is that the 
chem-milling process by which the sheets would be made has a +0. 0013 cm 
(±0.0005-in. ) thickness tolerance. This level of gradation produces local 
skin stiffness distortions of the order of + 15%. However, this figure can 
probably be improved upon, at additional cost, if several right and left 
cover- sheets are produced and the best two are selected for the model 

• Scaling of Standard Size Items - Consider the difficulty associated with 
reducing a standard size item such as a 1 3/4 inch bolt to l/lOth size. The 
nearest drill to produce a hole 1/I0th this size is a No. 16, which produces 
a 0. 450 cm (0. 177-in. ) hole with a 0. 010 cm (0. 004-in. ) tolerance. There- 
fore, either some level of compromise is required or higher costs to pre- 
serve the desired scale ratio are required 

• Assembly - Once the spar and rib members have been manufactured, the 
entire substructure, with the exception of the root and tip ribs, are positioned 
in a jig and spot welded together. Only one side of the skin corrugation is 
then welded in place. This is followed by welding the corresponding face- 
sheet using copper bars which slide in the spanwise direction through the 
corrugation crests, thereby providing back-up siqaports and electrical con- 
ductors for the welding operation. The second side of the fin is next assem- 
bled in the same order and manner. The assembly is completed by welding 
the root and tip ribs in place. Although machining of the spars, root joints, 
rudder actuator, and hinge pick-up fittings of the model appears feasible, 
many distortions from the prototype are expected, particularly those of the 
types presented in Figures 45-49 of Reference 8 in which small machining 
details were omitted or simplified 

3. 3 MODEL MANUFACTURING COSTS 

For a program involving a reduced scale replica model, manufacturing represents 
the major expense of the overall program. The reason for this is that the necessary 
miniaturization of members and gages requires elaborate handling procedures and 
attachment techniques requiring the development of special tooling fixtures, equip- 
ment handling, and fabrication methods. Table 3 summarizes the estimated 
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participation levels, by function, required for manufacturing the fin model. The 
material cost for the necessary amounts of rene 41 and inconel 718 are estimated at 
$6500, while the outside-vendor costs for chem-milling the variable-gage flat skins 
have been approximated at $50,000 for two complete sets. It should be noted that 
all labor time estimates have approximately a ±20% tolerance associated with them. 
Furthermore, it is estimated that an additional 5 man-months of engineering design 
and stress analysis would be required, beyond the level demonstrated by Fig. 10, to 
produce a complete set of drawings for the 1/I0th scale fin model. 
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SECTION 4 THERMO-STRUCTURAL ANALYSES AND MODELING ERRORS 


4.1 THERMAL ANALYSES 

4.1.1 Lumped Analysis of Entire Fin Structure 

The section of fin analyzed is bounded by the carbon/carbon leading edge, the 
aft-fuselage, and the rudder. These three interfaces are considered to be adiabatic 
surfaces in the thermal idealization. 

The fin thermal nodes were divided into cover, web, diagonal and stud-member nodes, 
as well as corner nodes, thus making them correspond almost identically to the structural 
model members and nodes (Fig. 8). The node locations for a typical chordwise segment are 
shown in Fig. 13. Prior to earth touchdown, a "space sink" at -273°C (-460°F) was rep- 
resented by a fictitious node. After touchdown, the external surfaces of the fin were con- 
vectively coupled to ambient air at 20°C (70°F) through a second fictitious node. 

Thermal properties and heat transfer coefficients used in calculating thermal 
capacitance, and conductive, convective, and radiative couplings are listed in 
Table 4. 

TABLE 4 THERMAL PROPERTIES AND COUPLINGS FOR THE HOT FIN NETWORK 


Material 

Thermal Capacity, c 

imic< Btu/lb -° F > 

Density, 

gm o 

S’ (Ib/in. 
cm ' ' 

Thermal 

Conductivity, k 
g-cal 

m-sec-°C 
(Btu/sec-in. -°F) 

Rene 

Inconel 

0. 12 (0. 12) 
0. 12 (0. 12) 

8.23 (0.297) 
8.23 (0.297) 

3.45(1.93 x 10~ 4 ) 
3.45(1.93 x 10- 4> ) 

rr— cal 

Heat transfer coefficient between Fin/Ambient Air, h= 0.367 — 

Hr-cm 2 -°C 

Panel emissivity, E=0.8 (. 75 Btu/Hr-ft^-°F) 
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-V-BAR THERMAL NODES 
• CORNER THERMAL NODES 
■ PANEL THERMAL NODES 



Note: For clarity, some nodes 
have been left unnumbered 
in the sketch 


Fig. 13 Thermal Structural Idealization of H3T Orbiter Fin Between Ribs 9 and 10 
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Preliminary analyses of the H3T configuration had shown that maximum temper- 
ature gradients would occur in the fin during reentry. Therefore, this condition was 
analyzed in detail. Reentry heat flux as a function of time for a 650°C (1200°F) 
isotherm is shown in Fig. 14. Heat flux levels at other locations of the fin were 
estimated using a fourth -power-law ratio. For example, the heat flux at the 540OC 
(1000°F)* isotherm, Q^ooO’ was determined by using the relation: 

q 1000 = q 1200 00 \ 

\ 1200 / 

where 

T 1000 = 810C>K ( 1460 ° R ) T i2oo = 92oOK (1660 o R) 

Using an initial uniform temperature of 71°C (160°F) and reentry heat flux 
boundary conditions and the thermal network described, thermal analyses were per- 
formed for both prototype and 1/lOth scale model. Detailed digital results for each 
60-second time interval for the prototype and 6 -second interval for the model were 
obtained, and temperature history plots were made for each prototype structural 
element between ribs 9 and 10 (see Fig. 8). Once the required scale shift factor of 
10 is made in prototype to model time, the prototype and model plots were found to 
overlap for the scale used in Fig. 16, 17, 18, and 19. The results for panels 61, 

63, 65 and 67 are presented in Table 5 and the associated scaling temperature errors 
are plotted in Fig. 20. 

Since these results are typical of other panels throughout the idealized fin, 
a more careful search of the digital printout was made to determine the greatest 
error in nodal temperature throughout the structure. Points A and A' of Fig. 8, 
which are located at the intersection of a rib and the front spar, indicated the 
greatest discrepancies. Temperature history plots for these points in the model 
and prototype are compared in Fig. 21. The results may be interpreted as follows: 
since only radiation was scaled accurately, the largest scaling errors occur during 
periods of maximum conduction to interval radiation flux ratios. Because this ratio is 
approximately inversely proportional to absolute temperature level to the third power, 
thermal modeling errors are lowest at times corresponding to the higher tempera- 
ture levels (e. g between 40 and 180 seconds of model time in Fig. 21). Another 

*See Fig. 15 for the 540°C (1000°F) isotherm location. 
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TEMPERATURE 



TIME. SEC 


Fig. 16 H3T Fin Thermal Analysis — Cover Nodes 61, 63, 65 and 67 



Fig. 17 H3T Fin Thermal Analysis — Spar Nodes 347, 348, 349 
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TEMPERATURE TEMPERATURE 


195 




0 800 1600 2400 3200 4000 

TIME. SEC 

Fig. 18 H3T Fin Thermal Analysis — Vertical Truss Nodes 192 Through 196 



Fig. 19 H3T Fin Thermal Analysis — Diagonal Truss Nodes 285 through 288 
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Fig. 20 Cover Panel Temperature Errors in Hot Fin Model 
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Fig. 21 H3T Hot Fin Thermal Analysis — Maximum Error Points A and A' 
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COVER SKIN TEMPERATURES - °C (°F) - IN PROTOTYPE 
AND MODEL OF H3T ORBITER FIN 
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significant fact is that the largest errors occur at nodes with the largest conduction 
couplings. Since the analysis is based upon a discrete idealization of properties 
which are actually distributed, it is felt that the above errors are conservative and 
represent upper limits to those anticipated in an actual hardware test. 

4.1.2 Detailed Analysis of Corrugated Covers 

Three main sources of temperature-error exist in the model’s corrugated 
covers. These are due to the (1) difference in model time-scale between radiation 
and conduction heat transfer of a reduced scale model, (2) distortions required by 
the weld-flats in the model, and (3) gage-size tolerances which have a much greater 
effect upon the model than the prototype. 

The effect of each of these was evaluated by conducting simulated reentry 
analyses for separate thermal idealizations (see Fig. 22) of a typical segment. 
Analysis of the following four thermal idealizations were performed: 

• Prototype covers 

• l/10th size perfect replica model of covers 

• 1/lOth size distorted model (corrugations flattened at crest and trough) 

• l/10th size distorted model with maximum material gage size error 
tolerance 

Results of the four analyses are presented in Table 6 and Fig. 23 through 25. 
They indicate that the temperature errors are slight, major errors are associated 
with earlier reentry times, and that they are attributable to the inevitable conduction/ 
radiation model time-scale discrepancy. 

4.2 STRUCTURAL ANALYSES 

4.2.1 Lumped Analysis of Entire Fin Structure 

Since the maximum temperature errors between replica model and prototype 
are of the order of 5%, application of the scaling laws indicates that thermal strains 
and stresses will be almost identical. Hence, they were only calculated for the pro- 
totype. 

Preliminary evaluation of the maximum temperature gradient between ribs 9 
and 10 from the curves of Fig. 16 through 19 indicated that it occurred within 600 
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TEMPERATURE 


PROTOTYPE & 10 x REPLICA MODEL 
10 x FLAT CORRUGATION MODEL 


©V 


Fig. 22 Hot Fin Cover Thermal Idealizations 


NOTE: NODE 1 IS 

DEFINED IN FIG. 22 


PROTOTYPE 
REPLICA MODEL 



PROTOTYPE 0 100 200 300 400 500 600 700 800 9 


MODEL 0 10 20 30 40 50 60 70 80 9 

TIME, SEC 

Fig. 23 Hot Fin Cover Thermal Analysis — Node 1 



TABLE 6 COMPARATIVE TEMPERATURE (°C(°F)) RESULTS OF COVER SKIN THERMAL 

ANALYSIS FOR HOT FIN REENTRY 
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*See Fig. 22 for identification of node number 
**Maximum gage size error. 

***A11 results after 480 sec are identical for each model 





PROTOTYPE 

REPLICA MODEL 


NOTE: NODE 3 IS 

DEFINED IN FIG. 22 


I I I I I I I I 1 1 1 J 
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TIME, SEC 

Fig. 24 Hot Fin Cover Thermal Analysis — Node 3 
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TIME , SEC 


Fig. 25 Hot Fin Cover Thermal Analysis — Node 7 


seconds after reentry. A more precise evaluation of the maximum gradient condition, 
using computer printout, revealed that the maximum chordwise gradient occurred at 
approximately 300 and 420 seconds after reentry. Therefore, finite element internal 
thermal load analyses were performed for the 300- and 420-second time-slice temp- 
erature distributions. Because of necessary assumptions associated with the finite 
element method the covers, which consist of outer flat sheets attached to corrugated 
inner sheets, were idealized as single orthotropic membrane elements. The results 
showed no significant differences in internal load intensities between 300 and 420 sec- 
onds. Furthermore, these loads were an order of magnitude lower than those caused 
by the peak mechanical air loads used for design purposes. Similar to the mechanical 
loading cases, the thermally induced primary stresses in the orbiter fin were in the 
spanwise direction, although secondary loads also occurred in the chordwise and 
thickness directions. Inasmuch as the material allowables at peak reentry temperatures 
are only reduced by approximately 20%, the thermal stresses do not appear critical 
in themselves. 

Fig. 26 shows the thermally induced average spanwise coverloads between ribs 
9 and 10 computed by the finite element program and classical beam theory. As may 
be seen, the loads determined by beam theory are quite accurate and slightly conser- 
vative. This is to be expected as a key assumption of beam theory is that sections 
normal to the spanwise axis remain plane while heated. Since this assumption intro- 
duces more severe constraints than actually exists, the more accurate finite-element 
results are slightly lower and permit the calculation of thickness and chordwise stres- 
ses. However, the relative simplicity and accuracy associated with the use of ele- 
mentary beam theory for the entire cross section justified its use for obtaining the 
stress results which follow. 

A comparison of spanwise loads in corrugated and noncorrugated spar web 
designs for the 300 second time slice are presented in Table 7. The results clearly 
demonstrate how the corrugations greatly releave the loads induced in the forward 
(hotter) external skins as well as in the spar webs themselves. 

4. 2. 2 Detailed Analysis of Corrugated Skin 

A more detailed thermal analysis of the cover (see Fig. 22) between ribs 9 
and 10, which permitted the determination of temperature differences in the built-up 
construction, revealed that the 300-second time slice produced more severe cover 
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TABLE 7 AVERAGE STRESSES IN CORRUGATED AND 
NON CORRUGATED H3T HOT FIN SPAR-WEB 
300-SECOND TIME SLICE 


Element 

Number 
(See Fig. 13) 

** 

Case 

Average 

Stress* ** 

newtons /cm^ 

(psi) 

Front 


C 

0 

(0) 

Spar Web 

347 

N 

65, 700 

(95, 250) 

Front 


C 

51, 100 

(74, 100) 

Spar Cap 

623 

N 

40,700 

(59, 000) 

Front 

61 

C 

-13,100 

(-19,000) 

Cover 


N 

-22, 100 

(-32,000) 

Second 

63 

C 

-4,490 

(-6,500) 

Cover 


N 

-11,000 

(-16, 000) 

Third 

65 

C 

6,900 

(10, 000) 

Cover 


N 

1, 380 

(2, 000) 

Mid -Spar 

348 

C 

0 

(0) 

Web 


N 

24,500 

(35,500) 

Mid -Spar 


C 

26,900 

(39,000) 

Cap 

629 

N 

22,000 

(32, 000) 

Aft 

67 

C 

3,450 

(5,000) 

Cover 


N 

0 

(0) 

Aft 


C 

0 

(0) 

Spar Web 

39 

N 

32,600 

(47,300) 

Aft 


C 

26,400 

(38, 250) 

Spar Cap 

631 

N 

24,400 

(35,400) 


*Does not include stress resulting from temperature variation through the 
thickness of the cover skin. These latter self-equilabrating stresses must 
be superimposed upon the above cover stresses. 

**C = corrugated 
**N = noncorrugated 
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stress than the 420 second time slice. Thus, thermal stresses in the covers between 
ribs 9 and 10 are presented only for the 300-second time condition (see Table 8). 

4.2.3 Effect of Heating Upon Bending and Twisting Stiffness 


Since the fin is free to expand in the spanwise direction, the effect of temperature 
upon the bending stiffness is primarily due to a reduction in material allowables which, 
for the materials and temperatures considered, is under 20%. 

An equation relating the change in torsional stiffness of a beam-type structure 
during heating, which causes a reduction in material properties and axial stresses, 
is (see References 9 and 10): 

<“>e£t <“>H Ja 17 x * 2 dA 

= + . 

<GJ ) 0 (GJ) q <GJ ) 0 

where GJ are twisting stiffnesses of the cross-sections, the a x are spanwise stresses, 
and the subscripts "o" and H correspond to unheated and heated conditions, res- 
pectively. 

The twisting stiffnesses were computed by modifying the thin walled multi-cell 
torsion formula given in Reference 11 to include a variable shear modulus, G, which 
depends upon temperature. Results for the corrugated shear web design are: 


(GJ) q = 117 x 10 3 newtons - m 4 (63.2 x 10 9 lb-in. 4 ) 

(GJ) h 

— - — x 100 = 92. 5% 

< GJ >o 

Sa ^x r2 

— x 100 = - . 25% 

GJ o 


(GJ) 


feff 


(GJ) 


x 100 = 92. 2% 


It should be observed that the effect of spanwise stress upon torsional compliance is 
extremely small in this case while the reduction in material properties causes the primary 
change in twisting stiffness. 
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TABLE 8 H3T FIN MID-PANEL SKIN STRESSES BETWEEN RIBS 9 AND 10 
(REF. FIG. 8 ) AT 300 SECONDS AFTER REENTRY 
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y Left 
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SECTION 5 THERMAL DYNAMIC TESTING 


A practical study into various aspects of testing the model described and 
analyzed in Sections 3 and 4 was performed. Special consideration and highlights of 
this effort, as they apply to the physical realization of the hot fin and thermal dynamic 
model testing in general, arc described below. 

5.1 METHODS AVAILABLE 

For geometrically reduced models made of the same material as the prototype, 
the discussion of Section 2. 1 indicates that the time scale, r , must shrink in exactly 
the same ratio as the linear dimension scale, X , for problems in which radiation is 
the dominant mode of internal heat transfer. Thus, to simulate fin healing associated 
with reentry, a 1/1 Oth scale model would have to achieve the same temperature levels 
in a time span of the order of three minutes. Thus, the structural properties will be 
changing rapidly. As possible solutions to the problems this introduces, a number 
of techniques for obtaining natural vibration frequencies and mode shapes during a 
short time span are presented. 

5.1.1 Single Mode Methods 

Two mode-by-modc procedures to obtain natural vibration frequencies and mode 
shapes during transient heating are presented in References 12 and 13. The methods 
discussed attempt to electronically track a single mode as its frequency varies (due 
to transient heating) by means of a servo control which constantly adjusts the input 
frequency to maintain the system in resonance. Depending upon the type of electronic 
circuitry employed, the control systems which have been used respond to either peak 
amplitude or phase angle. A third approach used in this family of single mode methods 
employs a fixed excitation frequency during heating, while the output instrumentation 
records the precise instant when the system is in resonance at the applied forcing fre- 
quency. Needless to say, these test procedures require many repeated applications 
of the thermal duty-cycle to track a reasonable number of structural modes, as well 
as the servo electronics to effect the desired control. 
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5.1.2 Multimode. Methods 

Alternate and potentially very attractive procedures for determining frequencies 
and mode shapes of structures with rapidly changing properties involve the processing 
of transient data to obtain steady-state information. The main differences among the 
various methods, as will be seen, is in the type of structural dynamic excitation re- 
quired. 

The "Impulse Response Method" is based upon measuring the response history, 
hjj(t), at various structural locations, i, to an impulse load applied at point j. By 
comparing the Fourier transformed equations of motion to impulsively loaded structures 
with those for sinusoidal excitation, it can be shown that a complex frequency response 
at point i, H„ ( oj ), caused by steady-state vibration excitation applied at j, may be 
calculated from the relationship 

r°° _ 

H ( W ) =j h (t) e -laJt dt 
00 J 


Thus, H_ may be evaluated numerically once the h^ are determined from recorded 
experimental data. For lightly damped structures with reasonably spaced frequencies, 
natural frequencies and mode shapes may be extracted from these calculations since 


where ^ 

A number of practical aspects associated with implementation of this method 
have been investigated and documented in a series of reports issuing from the Uni- 
versity of Southampton (References 14 through 16) and Grumman Aerospace (Ref- 
erences 17 through 19). Problems associated with excitation level and duration, 
length of recorded output data, closely spaced natural frequencies and data reduction 
using proper filtering and interpretation techniques, are discussed in these works 
with reference to tests actually performed. 


<t>. 

i 


<k) 


H ij k> 


| is the modal vector corresponding to the kth natural frequency, cu 


The testing procedure used in the "Rapid Frequency Sweep Method" is to load 
the structure by a constant amplitude excitation of rapidly varying frequency (Ref- 
erences 14 and 20). The idea behind this is similar to the impulse method in that the 
Fourier transforms of input and response may be used to compute the complex fre- 
quency response obtained from conventional steady-state vibration testing. 
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Advantages of the rapid frequency technique over the impulse response method 
are: the input energy can easily be made larger since a longer exicitation time is in- 
volved (hence, the output is also larger and less susceptible to measurement and pro- 
cessing error) and the range of resonances excited can be more accurately controlled, 
thus simplifying the required analogue-to-digital conversion of response data. As 
compared with steady -state testing techniques, transient tests, such as the impulse 
response and rapid frequency sweep methods minimize testing time, but greatly in- 
crease the data processing task. 

The "Step Input Method" is similar to the impulse and rapid frequency sweep 
test methods in that it uses Fourier transform computations of recorded transient 
structural responses, except that now the excitation input is a suddenly released static 
displacement of the structure. This technique was proposed by R. Mazet (Reference 
21) and was also used with success by others (References 22 and 23) 

5.2 ADDITIONAL TEST CONSIDERATIONS 

5.2.1 Heating Inputs 

The surface area involved and thermal flow rates required for providing reentry 
temperatures in the present study are sufficiently large so as to preclude the use of 
aerodynamic wind tunnel heating. This suggests the use of radiant heating by quartz 
lamps which are programmed and controlled to deliver desired temperature-time 
curves over various regions of the structure. Each region is designed to have an 
attached thermocuple follow a specified temperature history. The problem of heat 
leakage at boundaries can be overcome by oversized heating elements or by having separate 
adjacent regions controlled by guard heaters. 

Because of the requirement of high temperatures for short times, it is possible 
that the cooling phase of the curve may be the most difficult to simulate thermally 
For example, the hot fin skin temperature drops some 500C° (900F°) in 100 seconds (see 
Fig. 21). Thus, the test may well require the use of cool air blowers and baffles. In 
connection with this, it may be desirable to have more precise control elements. A 
laboratory facility which has a heating power capability of some 500 kva should be adequate 
for the combination of temperature, mass and thermal rates required for the hot fin. 
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5.2.2 Instrumentation 


The principal types of instruments available for thermal dynamic testing are 
thermocouples, accelerometers and strain gages. Besides being capable of funct- 
ioning at elevated temperatures, these instruments should be light-weight, accurate, 
and capable of withstanding the dynamic environment. 

The Ceramo Co. markets a 0.1cm (1/25 in.) OD Chromel-Alummel thermo- 
couple wire with an upper temperature limit of 760°C (1400°F). In addition, it is purported 
to have a fast response time and a ±3/4% error tolerance. Light-weight electro- 
formed or etched Ni-Cr strain gages which can be welded in place are sold by 
Microdot Inc. These are conservatively rated for operation up to 510°C (950°F) 
but have been reported for use in environments of 680°C (1250°F). In addition. 

Microdot also makes a platinum tungsten weldable alloy strain gage for dynamic 
measurements up to 815°C (1500°F). 

Small and light accelerometers (1.27cm (1/2 -in.) OD and 20 gm (0.7 oz), res- 
pectively) designed for use at temperatures upto760°C (1400°F) are produced by Endevco. 
However, they are expensive (approximately $1,000 each) and, because of their low 
gain at the temperatures involved, must be operated in conjuction with special charge 
amplifiers (a possible additional expense). Thus, they should be used more sparingly 
than strain gages which are much less costly, can be used to obtain comparable 
natural frequency data and, with some additional effort, may be calibrated to give 
modal data. 

A relatively new method for measuring vibtation of heated structures employs 
pulsed differential holographic measurements to produce modal interference fringe 
patterns (References 24 through 26). Although initial acqusition of such a system is 
costly (somewhere between $50 and $100K) it could provide useful for structural 
temperature levels in excess of 1100°C (2000°F). In addition, the method is non- 
contacting and produces entire modal patterns rather than data for a limited number 
of discrete points. 

5.2.3 Mechanical Loading 

Long rods or reciprocating arms can be used to transmit loads to specimens in 
the heating area. These can be connected directly to shakers (Reference 12) and to cam 
mechanisms designed to deliver an impulse or sudden load-release (References 10 and 
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22). In many cases, attention must be given to ensure that serious thermal or 
mechanical distortions are not introduced at the excitation points. 

If stiffness measurements during heating are required, or if programmed dynamic 
loads must be applied simultaneously in the heated area, piezo-electric quartz rods 
can be used in conjunction with servo control and potentiometers to regulate the 
mechanical loading with an electric input. 

5.2.4 Thermal Cycling 

The effects of thermal cycling associated with repeated testing of the specimen 
should be considered since inelastic deformations might have the effect of changing 
the structure from test to test. This is most likely to happen in the first few thermal 
cycles as specimen heating tends to release the manufacturing stresses through de- 
formation associated with material annealing. Since the structures are ordinarly 
under a 1-g loading only and primary structure materials are normally rather stable, 
it is expected that the importance of creep phenomena will be small. 

5.2.5 Cost Estimates For Testing the Hot Fin 

Table 9 lists the man-months and type engineer or technician responsible for 
each experimental task associated with the hot fin test (cost estimates for model 
design and manufacture are presented in Section 3. 3). Instrumentation charges are 
additional and are estimated at $2,000 for two high-temperature accelerometers 
and another $2, 000 for thermocouples and high temperature weldable strain gates. 

Data reduction, involving conversion of analogue tapes to digital data, will require 
approximately another two man-months of effort. 

The experimental program estimates are based upon the following set of exper- 
iments: thermal tests without mechanical excitation, dynamic tests at room temp- 
erature, and combined thermal-dynamic tests. The unheated dynamic tests were 
assumed to consist of a frequency sweep and one transient type (i. e. sudden step or 
impulsive loading test). The purpose of the room temperature dynamic tests will 
be to provide data for assessment of the thermal dynamic test data through a com- 
parison of transient and steady-state results. 
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ESTIMATED MAN-MONTHS FOR TEST OF l/lOth SCALE HOT FIN MODEL 
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SECTION 6 CONCLUSIONS AND RECOMMENDATIONS 


A search of the literature dealing with thermal effects upon dynamic response of 
engineering structures has uncovered a number of references concerned primarily with 
either analysis (References 9, 27 and 28) or test (References 10, 22 and 23), but none 
which correlate the two in a well coordinated and detailed manner. Another observation 
made, and concurred with in a recent note by Professor Swaim (Reference 29), is that 
this technical area of investigation was active in the late fifties and early sixties, but very 
little has since been reported upon. Consequently, there appears to be a need in the 
shuttle effort to perform concurrent thermal-structural tests and analysis programs and 
to evolve appropriate confidence levels to be used in evaluating theoretical predictions. 
Full-scale tests should be performed upon prototype details, such as corrugated skin- 
stringer attachments and mechanical joints, but larger components such as wings and 
tails will require the use of reduced scale models. In connection with such experiments, 
the scaling laws of Section 2 and the instrumentation and test-method information of 
Section 5 may be of assistance. 

Although tests can be performed at room temperature upon models with manu- 
factured reduced stiffness, thus simplifying the need for simultaneously supplying heat 
and mechanical excitation, such models would be severely restricted, corresponding to 
only one time instant of a single maneuver. Thus, transient heating tests which simulate 
an entire environmental phase are preferred. 

In connection with such testing, it is felt that thermal scaling errors can be kept 
small but the possibility of shifting materials and temperature levels is quite impractical 
because of the number of important variables involved which are temperature dependent. 
Furthermore, the primary-structure temperature levels are not sufficiently high nor are 
the dynamic test techniques so difficult that they cannot be adequately managed in an 
acceptable modeling program with state-of-the-art instrumentation. 

Because of the high level of complexity and cost associated with reduced-scale 
replica models (witness the results of the manufacturing cost study upon the hot fin 
presented in Section 3) "adequate" thermal-dynamic models, rather than true-replica 
reduced scaled models, are recommended where possible. Of course this tack requires 
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an a priori knowledge of the basic parameters, since it assumes that unanticipated 
problems will not occur. If such phenomena did occur, they could adversely affect the 
overall stiffness modeling of a non-replica component in a way which did not correlate 
with behavior of the prototype. For example, if the linear bending to twisting stiffness 
ratios were identical in model and prototype, i.e. , 

El (El)' 

GJ (GJ)' 

but thermal stresses caused an unanticipated buckle in a reinforced cover skin (in 
either model or prototype but not both), the true stiffness ratios would not agree. 

A replica model would not have this deficiency since nonlinear phenomena which are 
strain or stress dependent (e.g. , buckling), but not time dependent (e. g. creep), 
would scale one to one. 
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APPENDIX 


THERMAL DYNAMIC SCALING LAWS FOR A BEAM 


The governing equation for transverse beam vibrations under the influence of 
heating and axial loading is 


dx 
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El -M t 

+ — 
d x 

(EA € x -N t ) 

a> a> 
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where w is the transverse beam deflection, z and x are the transverse and axial beam 
coordinates, E, A and p are the beams elastic modulus, area and mass density, 
respectively, and t is time. The remaining variables are defined below. 
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stretching stiffness 
mass per unit length 

strain displacement relationship 


where u is the centroid 's axial deflection, € is the axial strain and a is the coef- 

X 

ficient of thermal expansion. 

An additional assumption required in the derivation of the vibration equation in- 
volves definition of the beam centroid (at each temperature state) which is defined by 

/ez dA = 0 

Introducing the scaling parameters, 
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where the primed and unprimed, quantities relate to the model and prototype, respectively Substitu- 
ting the above equations into tlif nonlinear strain-displacement relationship for the model, yields 
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Recognizing that it is desirable to have the model strain, e x , proportional to 


the prototype strain, yields 
X u 


x 

Thus, 


€ x = and t_he remaining auxiliary scaling laws become 
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The main scaling laws, reduce to 




and 
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where 

It should be noted that axial vibration effects have been ignored in the above de- 
rivation. If they were included, the last set of equations would result in 


X x =X , 


a 6 = a and T 
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